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A  transient  two-phase  mass  transport  model  for  liquid  feed  direct  methanol  fuel  cells  (DMFCs)  is  devel¬ 
oped.  With  this  model,  various  processes  that  affect  the  DMFC  transient  behaviors  are  numerically  studied. 
The  results  show  that  the  cell  voltage  exhibits  an  overshoot  behavior  in  response  to  a  sudden  change  in 
the  current  density.  The  magnitude  of  the  overshoot  depends  on  the  magnitudes  of  the  change  in  the 
cell  current  density  and  the  initial  current  density.  It  is  found  that  the  dynamic  change  in  the  methanol 
permeation  through  the  membrane  to  the  cathode  results  in  a  strong  cathode  overpotential  overshoot, 
which  is  believed  to  be  the  predominant  factor  that  leads  to  the  cell  voltage  overshoot.  In  contrast,  the 
anode  overpotential  is  relatively  insensitive  to  the  changes  in  the  methanol  concentration  as  well  as  CO 
surface  coverage  in  the  anode  catalyst  layer.  Moreover,  the  effect  of  the  double  layer  capacitance  (DLC) 
on  the  cell  dynamic  behavior  is  studied  and  the  results  show  that  the  DLC  can  smoothen  the  change  in 
the  cell  voltage  in  response  to  a  change  in  the  cell  current  density.  Furthermore,  the  dynamic  response  of 
mass  transport  to  a  change  in  the  cell  current  density  is  found  to  be  rather  slow.  In  particular,  it  is  shown 
that  the  slow  response  in  the  mass  transport  of  methanol  is  one  of  the  key  factors  that  influence  the  cell 
dynamic  operation. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Because  of  its  promising  virtues  for  powering  portable  electronic 
devices,  the  direct  methanol  fuel  cell  (DMFC)  has  been  extensively 
studied  over  the  past  decade,  most  of  them  focused  on  the  fuel  cell 
steady-state  behaviors.  As  a  matter  of  fact,  the  transient  behaviors 
of  the  fuel  cell  associated  with  the  change  in  operating  conditions, 
such  as  system  start-up,  shut-down,  sudden  changes  in  the  power 
level  are  also  important  and  need  to  be  investigated.  A  change  in 
the  cell  operation  conditions  may  result  in  the  changes  in  the  cell 
temperature  and  the  reactant  concentrations  at  the  electrode  sur¬ 
face,  which  eventually  determine  the  overall  transient  response  of 
the  cell.  A  better  understanding  of  the  fuel  cell  dynamic  charac¬ 
teristics  under  the  transient  operation  conditions  is  thus  essential, 
which  enables  the  engineers  to  better  design  and  control  the  fuel 
cell  system.  With  this  purpose,  some  experimental  investigations 
of  the  DMFC  dynamic  operating  characteristics  have  been  reported. 
A  widely  used  experimental  method  is  to  make  a  step  change  in 
the  input  parameter  to  the  DMFC,  such  as  the  methanol  feed  con¬ 
centration  or  the  current  density,  and  to  study  the  response  of  the 
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cell  voltage.  For  instance,  Sundmacher  and  co-workers  studied  the 
dynamic  response  of  a  DMFC  with  a  step  change  in  the  methanol 
feed  concentration  [1,2].  Other  researchers  also  investigated  the 
dynamic  behaviors  of  the  DMFC  by  making  a  step  change  in  the  cell 
current  density  [3-5].  In  response  to  the  change  in  the  cell  current 
density,  the  cell  voltage  usually  exhibited  an  overshoot  behavior. 
Argyropoulos  et  al.  [3,4]  also  studied  the  DMFCs  response  to  the 
consecutive  change  in  the  current  density  in  addition  to  the  single 
step  change  in  the  current  density.  In  addition,  Kallo  et  al.  [5]  stud¬ 
ied  the  cell  voltage  response  of  a  gas  feed  DMFC  to  a  step  change 
in  the  cell  current  density.  They  also  investigated  the  effects  of  the 
double  layer  capacitance  (DLC),  the  formation/removal  of  CO  poi¬ 
son  on  Pt  catalyst  at  the  anode,  and  the  methanol  crossover  on 
the  cell  voltage  response.  Besides  above-mentioned  works,  in-situ 
flow  visualization  studies  [6]  were  also  conducted  to  investigate  the 
transient  behaviors  of  two-phase  flow  in  the  flow  field  of  the  DMFC 
and  its  effect  on  the  cell  performance.  It  should  be  mentioned  that 
electrochemical  impedance  spectroscopy  (EIS)  can  also  be  used  to 
study  the  dynamic  cell  behavior.  Although  EIS  can  cover  a  wide 
time  range  as  well  as  monitoring  of  most  processes  occurring  in 
the  fuel  cell,  the  quantitative  interpretation  and  EIS  modeling  of 
the  complete  DMFC  have  not  been  achieved  yet. 

In  addition  to  experimental  efforts,  mathematical  modeling  can 
be  very  useful  in  understanding  both  the  steady-state  and  the 
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Nomenclature 

Alg  interfacial  specific  area  between  liquid  and  gas 
phase  (m2  m-3) 

As  specific  surface  area  of  the  active  reaction  sites 

(m2  m-3) 

C  molar  concentration  (mol  m-3 ) 

CAcl  double  layer  capacitance  of  anode  catalyst  layer 

(Cm-2) 

CCcl  double  layer  capacitance  of  anode  catalyst  layer 

(Cm-2) 

D  diffusivity  (m2  s_1 ) 

F  Faraday  constant  (96,478  C  mol-1 ) 

i+  proton  current  vector  in  the  catalyst  layer  (A  nrr2 ) 

1  proton  current  vector  in  the  membrane  (A  m-2 ) 

/ceii  cell  current  density  (A  nrr2 ) 

/para  parasitic  current  resulting  from  methanol  crossover 

(Am-2) 

jo  exchange  current  density  (A  m-2 ) 

ja  anode  current  density  (A  m-3 ) 

jc  cathode  current  density  (A  nrr3 ) 

/<ai  reaction  rate  constant  for  anode  reaction  step  1 
(molm-2  s_1) 

/<a2  reaction  rate  constant  for  anode  reaction  step  2 
(molm-2  s_1) 

kc  condensation  rate  (mol  atm-1  s_1  m-3 ) 

ke  evaporation  rate  (1  atm-1  s_1 ) 

/<h  Henry’s  law  constant  (Pa) 

kr  relative  permeability 

I<  permeability  of  porous  material  (m2 ) 

m  source  term  in  mass  conservation  equation 

(kgnrr3  s_1) 

M  molecular  weight  (kg  mol-1 ) 

nd  electro-osmotic  drag  coefficient 

N  vector  flux  of  species  (mol  m-2  s_1) 

pc  capillary  pressure  (Pa) 

pg  gas  phase  pressure  (Pa) 

Pi  liquid  phase  pressure  (Pa) 

R  gas  constant  (J  mol-1  K_1 ) 

R  source  term  in  species  conservation  equation 

(molm-3  s_1) 

R  interfacial  species  transfer  rate  (mol  m-3  s_1 ) 

/?!  reaction  rate  of  anode  reaction  step  1  (mol  m-2  s_1 ) 

R2  reaction  rate  of  anode  reaction  step  2  (mol  m-2  s_1 ) 

^contact  Ohmic  contact  resistance  ( £2m2 ) 
s  liquid  saturation 

T  temperature  (K) 

ug  velocity  vector  of  gas  phase  (m  s_1 ) 

Ui  velocity  vector  of  liquid  phase  (m  s_1 ) 

V0  thermodynamic  equilibrium  voltage  (V) 

Veen  cell  voltage  (V) 

x  coordinate,  m,  or  mole  fraction  in  liquid  solution 

(mol  mol-1) 

y  coordinate,  m,  or  mole  fraction  in  gas  mixture 

(mol  mol-1) 

Greek  letters 

r  surface  concentration  of  Pt  (mol  m-2 ) 

a  gas-void  fraction  ( 1  -  s) 

aa  anode  transfer  coefficient  at  anode 
otc  cathode  transfer  coefficient  at  cathode 


8  thickness  of  porous  layer  (m) 

s  porosity  of  porous  medium 

r\  overpotential  (V) 

k  ionic  conductivity  of  membrane  (£2_1  m_1 ) 

fi  viscosity  (kg  m_1  s_1 ) 

6C  contact  angle  (°) 

Oco  surface  coverage  of  CO  species  on  the  anode  Pt  cat¬ 

alysts 

p  density  (kg  m-3) 

a  interfacial  tension  (N  m_1 ) 

Superscripts 

eff  effective  value 

in  inlet  condition 

ref  reference  value 

sat  saturated  value 

Subscripts 
A  anode 

ACL  anode  catalyst  layer 

ADL  anode  diffusion  layer 

C  cathode 

CCL  cathode  catalyst  layer 

CDL  cathode  diffusion  layer 

G  gas  phase 

L  liquid  phase 

Mem  membrane 

M  methanol 

MV  methanol  vapor 

W  water 

WV  water  vapor 


dynamic  behaviors  of  the  DMFC.  Most  of  the  existing  DMFC  models 
are  restricted  to  steady-state  condition  [7-20],  and  the  transient 
models  for  DMFCs  are  scarce.  Sundmacher  et  al.  [21]  developed  a 
model  to  study  the  voltage  response  to  a  change  in  the  methanol 
feed  concentration.  The  results  revealed  that  the  voltage  under¬ 
went  an  overshoot  in  response  to  a  step  change  in  the  methanol 
feed  concentration.  Schultz  et  al.  [22]  also  studied  the  voltage 
response  to  a  step  change  in  the  methanol  feed  concentration. 
Their  one-dimensional,  multi-component  mass  transport  model 
was  developed  based  on  the  Maxwell-Stefan  equation.  The  model- 
based  analysis  revealed  that  the  cell  voltage  overshoot  was  resulted 
from  the  different  time  responses  of  the  anode  and  cathode  overpo¬ 
tentials  to  the  change  in  the  methanol  concentration.  In  addition  to 
the  studies  with  a  step  change  in  the  methanol  feed  concentration, 
other  researchers  also  developed  mathematical  models  to  study 
the  dynamic  voltage  response  of  the  DMFC  to  a  change  in  the  cell 
current  density.  For  instance,  Krewer  and  Sundmacher  [23]  numer¬ 
ically  revealed  that  the  methanol  crossover  is  the  dominating  factor 
that  causes  the  cathode  overpotential  overshoot  and  then  the  cell 
voltage  overshoot.  This  result  is  consistent  with  the  experimental 
study  by  Kallo  et  al.  [5].  More  recently,  Krewer  et  al.  [24]  presented 
another  model  to  investigate  the  voltage  response  to  a  change  in 
the  cell  current  density,  and  showed  that  the  anode  reaction  mech¬ 
anism  was  the  main  physicochemical  phenomena  that  cause  the 
anode  overpotential  overshoot,  hence  the  cell  voltage  overshoot. 
However,  this  result  was  somewhat  contradictory  to  their  earlier 
study  [23].  It  is  also  worth  mentioning  that  a  few  transient  models 
for  studying  the  start-up  process  of  passive  DMFCs  have  also  been 
reported  [25,26]. 
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Fig.  1.  Schematic  of  the  model  domain.  Roman  numerals  refer  to  the  boundaries  of 
the  modeling  geometries. 

Our  literature  review  indicates  that  most  previous  transient 
transport  models  for  DMFCs  [23,24,27]  assume  that  the  transport 
of  methanol  through  the  anode  diffusion  layer,  catalyst  layer  and 
the  membrane  is  fast  enough  and  can  be  treated  as  a  quasi-steady 
process.  This  assumption  is  questionable  as  the  mass  transport  of 
species  is  usually  much  slower  than  the  charge  processes  and  elec¬ 
trochemical  processes  at  the  electrode  [27].  Also,  we  found  that  the 
interpretations  to  the  voltage  overshoot  behavior  by  different  mod¬ 
els  are  contradictory.  Motivated  by  addressing  these  problems,  in 
this  work  we  developed  a  transient,  two-dimensional,  two-phase 
mass  transport  model  for  DMFCs,  in  which  the  assumption  of 
quasi-steady  mass  transport  is  eliminated.  With  this  model,  var¬ 
ious  processes  that  affect  the  DMFC  transient  operating  behaviors 
are  numerically  studied. 

2.  Formulation 

Consider  a  two-dimensional  physical  domain,  as  sketched  in 
Fig.  1,  which  represents  the  typical  geometry  of  a  membrane  elec¬ 
trode  assembly  (MEA)  that  consists  of  anode  diffusion  layer  (ADL), 
an  anode  catalyst  layer  (ACL),  a  membrane  (MEM),  a  cathode  cat¬ 
alyst  layer  (CCL),  and  a  cathode  diffusion  layer  (CDL).  The  MEA  is 
sandwiched  between  two  parallel  flow  field  plates.  Since  both  chan¬ 
nel  width  and  rib  width  in  the  parallel  flow-field  are  symmetrical 
with  respect  to  their  middle  points,  only  a  half-rib  width  and  a 
half-channel  width  need  to  be  considered  to  save  the  computing 
cost. 

2.1.  Anode  porous  region 


Pi, a  :  j|(/Ol£Si,a)  =  v-  (piKjf  ■  V'Pl,a)  +  mu  (4) 

where  NM1  and  NM)g  are  both  vectors,  representing  the  trans¬ 
port  flux  of  methanol  and  methanol  vapor  in  the  anode  electrode, 
respectively.  Taking  into  account  the  main  transport  mechanism 
of  methanol,  namely,  the  molecular  diffusion,  macroscopic  con¬ 
vection  and  electro-osmotic  drag,  the  overall  flux  of  methanol  and 
methanol  vapor  can  be,  respectively,  calculated  from 

0/1,1**!  -  ^Mfl V  '  CM,1  ADL 

’  i  (5) 

CM,i«i-D^f,V-CMjl  +  ndjWxM-i  ACL 

and 

NM,g  —  Qvi,gUg  —  D^g  V  •  CM,g  (6) 

To  account  for  the  effect  of  methanol  evaporation  and  conden¬ 
sation,  the  interfacial  transfer  rate  of  methanol  between  the  liquid 
and  gas  phase  is  given  by: 

Rm  =  Alghlgs(  1  -  s)(P^r/Mv)  (7) 

where  denotes  the  saturation  pressure  of  methanol  vapor.  Note 

that  the  interfacial  transfer  of  methanol  between  the  phases  is 
finally  embodied  in  the  general  source  terms  on  the  right-hand 
sides  of  Eqs.  (l)-(4). 

2.2.  Cathode  porous  region 

The  two-phase  mass  transport  behaviors  in  the  cathode  porous 
region  (the  cathode  diffusion  and  catalyst  layer)  are  related  to  four 
major  variables,  i.e.,  oxygen  concentration  (C02  ,g ),  water  vapor  con¬ 
centration  (Cwv.g).  liquid  saturation  (s\)  and  gas  pressure  (pg).  The 
governing  equations  of  conservation  corresponding  to  the  four  vari¬ 
ables  are,  respectively,  given  by: 


c°2,g:  ^[£(1 -Si,c)Co2,g]  =  - 

-v-No2,g  +  *o2,g 

(8) 

9 

Cwv,g  •  -  sl,cKwV,g]  = 

:  -V  •  Nwv,g  +^o2,g 

(9) 

d  ,  .  _ _ _  r  _  _  kri 

Si,c  :  -gi(A£Si,a)=V  ■  P\l<-- 
dt  L  JL\ 

Vsl,c+V  Pg,c)  +dli5c 

(10) 

9  / 

Pg,c  :  g^[pg£(l -S,>c)]  =  V  -  ( 

•  V  Pg,cJ  +  dig?c 

(ID 

where  No2,g  and  Nwv,g  are  both  vectors,  denoting  the  transport  flux 
of  gas  oxygen  and  water  vapor  in  the  cathode  porous  region,  which 
can  be  obtained  from: 


In  the  anode  porous  region  (the  anode  diffusion  and  catalyst 
layers),  we  are  interested  in  four  major  variables,  including  liq¬ 
uid  methanol  concentration  (CM1),  methanol  vapor  concentration 
(Qvi.g),  liquid  saturation  and  liquid  pressure  (pi).  The  general  gov¬ 
erning  equations  corresponding  to  each  of  these  variables  are  given 
by 


Qvi,i 


.  d_ 
:  di 


(^i?aCM)i)  =  -V  •  NM)1  +  RM\ 


(1) 


Nj,g  =  Cj  gUg  —  Df gV  •  Cjjg  i:02,WV  (12) 

To  consider  the  effect  of  water  evaporation  and  condensation, 
the  interfacial  transfer  rate  of  water  between  the  liquid  and  gas 
phase  is  given  by  [28,29]: 


ke  (ywvPg  -  P W  )  y  WVPg  < 

kc  e(!  -  ^)ywv  (y  wvPg  _  psat^ )  ywvPg  >  psat^ 


(13) 


CM,g  =  g^[e(l  _sl,a)Qvi,g]  =  -V  •  NM,g  +^M,g 


(2) 


Si, a  :  7^[Pge(l-Sl,a)]=V. 


(tv  s'-+v ") 


'  dp, 


+dig,a  (3) 


where  p^  and  j/wv  denote  the  saturation  pressure  of  water  vapor 
and  mole  fraction  of  water  vapor  in  the  cathode  gas  phase,  respec¬ 
tively.  Note  that  the  interfacial  transfer  of  water  between  the  phases 
is  embodied  in  the  source  terms  on  the  right-hand  sides  of  Eqs. 
(8)— ( 11 ). 


1134 


W.W.  Yang,  T.S.  Zhao  /  Journal  of  Power  Sources  185  (2008)  1131-1140 


Table  1 

Constitutive  relations 


Parameters 
Capillary  pressure 

Relative  permeabilities 
Effective  diffusion  coefficients  of  species  [18] 


General  generation  rate  of  mass  in  liquid  phase 


General  generation  rate  of  mass  in  gas  phase 


Mole  generation  rate  of  species 


Expressions 


Pc  =  Pg  -  P\  =  cr  cos  6c(e/I<rbJ(s) 

f  1.417(1  -s)-  2.120(1  -s)2  +  1.263(1  - 
\  1.417s  -  2.120s2  +  1.263s3 
Liquid 


J(s)  = 


0<0C<  90° 

90°  <  0C  <  180° 


Kl=s 3 

kTg  =  (1  - 


Den  _  £)  c 

i.g  1,g 


5(1  -s)1 


Rq2  = 


dm,i8' 


Derr  =  . 


™l,a  =  • 


mu  =  • 


i :  02,WV,  MV 
(£  +  £N) 


ADL 

ACL 


[£/(DMje1-5S1-5)  +  Sn/Pm.nSn1'5)] 

Dm^nS1*5  MEM 

-MH2o^w  -  MmRm  ADL 

-(^h2o  +  ^  -  Mh2o^w  -  MmRm  ACL 

-Mh2o^w  CDL 

CCL 

1^h20^w  -  MmRm  ADL 

MCo2  -  MH2o^w  -  MmRm  ACL 


f  MH2o^w 


CDL 

-  MH2o^w  CCL 


1 

^  =  {  Sm  ■  «WV,a={? 


H 

4F 

-Rm 


RWV,c  = 


_  f  Rw  CDL 
“  \  Rw  CCL 


kw  ADL 
Rw  ACL 


2.3.  Membrane 


2.4.  Boundary  conditions 


Unlike  in  the  anode  and  cathode  electrode,  in  the  membrane, 
only  transport  of  liquid  phases  (dissolved  methanol  and  liquid 
water)  needs  to  be  considered  as  the  membrane  is  usually  regarded 
as  a  gas  insulator.  The  mass  conservation  of  dissolved  methanol  is 
described  by: 

9 

^(SmemCMeOH)  =  -V  •  NM  (14) 


As  indicated  in  Fig.  1,  the  computational  domain  is  enclosed 
by  six  boundaries.  The  boundary  conditions  at  each  boundary  are 
described  below. 

Boundary  t :  This  boundary  represents  the  inlet  of  reactant  sup¬ 
ply  at  the  anode,  at  which  the  concentration  of  liquid  methanol, 
the  concentrations  of  methanol  vapor,  liquid-phase  pressure  and 
liquid  saturation  are  all  specified  to  be  inlet  conditions: 


with  Nm  denoting  the  vector  flux  of  methanol  in  the  membrane. 
Generally,  transport  of  methanol  through  the  membrane  depends 
on  molecular  diffusion,  electro-osmotic  drag  and  convection,  and 
the  last  one  is  usually  negligible.  Accordingly,  the  flux  of  methanol 
in  the  membrane  is  given  by 

Nm  =  -^m[n  '  VCMeOH  +  nd,Mp  (15) 

With  respect  to  the  water  transport  through  the  membrane,  the 
flux  due  to  the  molecular  diffusion  can  be  ignored  if  the  membrane 
is  in  equilibrium  with  liquid  water  on  the  both  sides.  Hence,  the 
flux  of  water  across  the  membrane  (Nw)  is  given  by: 


NW  =  Hd,H2Op 


P\  K  Api  c_a 
Mh2o  Mi  $mem 


(16) 


Up  to  this  point,  we  have  presented  all  the  governing  equations 
that  describe  the  dynamic  two-phase  mass  transport  processes 
in  a  DMFC.  To  make  the  above  governing  equations  closed,  some 
constitutive  correlations  and  definitions  are  needed.  These  include 
capillary  pressure,  relative  permeability  for  both  gas  and  liquid 
phase,  effective  diffusion  coefficients  for  each  species  and  the 
source  terms.  All  these  correlations  and  associated  nomenclatures 
are  clearly  listed  in  Table  1. 


CM,1  =  CiJ,  CMV  =  CiJv,  s  =  l,Pl=pjn  (17) 

Boundary  2:  This  boundary  is  the  interface  between  anode  diffu¬ 
sion  layer  and  the  anode  rib  collector,  which  is  an  impermeable 
wall.  Accordingly,  all  the  fluxes  in  the  x  direction  are  zero  at  this 
boundary: 

^=0,  ^Y=o,  J=0,  fl=0  (18) 

OX  ox  ox  ox 

Boundaries  3  and  6:  These  two  boundaries  are  symmetrical,  respec¬ 
tively,  with  respect  to  the  middle  point  of  rib  width  and  the  middle 
point  of  channel  width.  Hence,  the  gradients  of  all  the  variables  in 
y  direction  are  set  to  zero: 

^  =  0  0  :  CM)1,  CMv,  s,  Pi,  Q)2,  Cwv,  (I9) 

Boundary  4:  Similar  to  Boundary  2,  Boundary  4  represents  the 
interface  between  cathode  diffusion  layer  and  the  cathode  rib  col¬ 
lector,  at  which  all  the  fluxes  in  the  x  direction  are  zero: 


9Q)2  _  n  dCwv 

dx  ~  ’  dx 


ds  9p2 
=  0,  =0,^=0 
OX  OX 


(20) 


Boundary  5:  This  boundary  represents  the  inlet  of  oxygen  supply 
and  outlet  of  water  removal  at  the  cathode,  at  which  the  following 
boundary  conditions  are  specified: 


Co2  =  C“  ,  Cwv  =  C^,s  =  0, Pg  =  Pgn 


(21) 
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The  conditions  at  the  interfaces  between  those  functional  lay¬ 
ers  in  the  ‘sandwiched’  cell  are  given  based  on  the  principle  that 
the  continuity  and  mass/species  flux  balance  are  required  at  each 
interface  to  satisfy  the  general  mass  and  species  conservation  of  the 
entire  cell.  The  more  details  can  be  found  elsewhere  [18]. 

2.5.  Electrochemical  kinetics 

On  the  DMFC  anode,  a  simplified  two-step  reaction  mechanism 
is  used  to  describe  the  electrochemical  oxidation  of  methanol  at 
the  platinum-ruthenium  catalysts  in  the  anode  [1].  In  the  first 
irreversible  step,  methanol  is  adsorbed  electrochemically  to  the 
platinum  catalyst,  forming  strongly  adsorbed  carbon  monoxide, 
four  protons  and  four  electrons. 

CH3OH  +  Pt^  CH3OHads-^COads+4H++4e-  (1) 

The  second  step  is  a  combination  of  two  consecutive  reactions: 
(a)  water  adsorbs  electrochemically  to  the  ruthenium  catalyst  in 
a  reversible  reaction  to  form  OHads,  a  proton  and  an  electron,  (b) 
the  formed  OHads  the  reacts  irreversibly  with  COads  to  form  carbon 
dioxide,  proton  and  electron.  In  these  two  consecutive  reactions, 
the  reaction  (b)  is  assumed  to  be  the  rate-determining  step,  while 
the  reaction  (a)  is  assumed  always  in  quasi-equilibrium.  Simply,  the 
overall  reaction  of  the  second  step  is  described  by: 

COads  +  H20R^2C02  +  2H+  +  2e“  +  Pt  (2) 

Since  both  the  reactions  steps  are  irreversible  electrochemical 
reactions,  their  rate  expressions  are  formulated  in  the  form  of  Tafel 
kinetics.  To  consider  the  effect  of  the  adsorption  and  desorption  of 
species  to  the  platinum  catalyst  surface,  the  rate  expressions  also 
contain  a  term  to  account  for  the  influence  of  the  surface  cover¬ 
age  of  the  platinum  catalyst  with  adsorbed  CO.  The  details  are  as 
follows. 

R]  -  fcal  ^f°H  (1  -  ^co)  exp  ( Q^F  >7a)  (22) 

CMeOH  V  RT  ’ 


difference  between  the  real  potential  and  the  thermodynamic- 
equilibrium  potential  of  ORR  which  are  also  referenced  to  a  SHE. 

In  addition  to  the  ORR,  a  second  undesired  reaction  also  takes 
place  in  the  CCL,  i.e.,  methanol  crossover  to  the  cathode  is  electro- 
oxidized.  This  parasitic  reaction  is  usually  very  fast  so  that  the 
methanol  crossover  is  immediately  and  entirely  consumed  in  the 
cathode.  Notice  that  this  parasitic  reaction  ultimately  results  in  the 
cathode  mixed  overpotential. 


2.6.  Cell  current  density  and  cell  voltage 


The  equivalent  circuit  for  a  DMFC  can  usually  be  simpli¬ 
fied  to  a  series  connection  of  the  membrane  resistance  and  two 
capacitor/resistance  parallel  circuits  corresponding  to  the  dou¬ 
ble  layer  capacitance  and  the  charge  transfer  resistance  for  the 
anode  and  cathode,  respectively.  Under  dynamic  operations,  the 
charging/discharging  process  of  the  DLC  may  occur,  accompa¬ 
nying  the  dynamic  variation  in  the  electrode  potential.  If  true, 
part  of  the  current  will  be  used  to  charge  the  DLC  or  released 
from  the  DLC.  Consequently,  the  dynamic  cell  current  density 
satisfies 


drja  _  .  fhcJ*dxdy 

dt  cel1  (Wc  +  Wr)/2 


(27) 


where  CAcl  represents  the  DLC  of  the  anode  electrode.  Note  the  first 
term  on  the  left-hand  side  of  the  Eq.  (27)  represents  the  current 
required  for  charging/discharging  of  the  DLC,  which  will  vanish  at 
the  steady-state  operation. 

As  mentioned  earlier,  methanol  may  permeate  through  the 
membrane  to  the  cathode.  Usually,  the  ‘parasitic’  current  density 
is  used  to  represent  the  rate  of  methanol  crossover  because  this 
part  of  current  cannot  be  used  to  do  work.  On  the  DMFC  cathode,  it 
is  assumed  that  both  the  cell  current  and  the  ‘parasitic’  current  are 
entirely  consumed  by  the  ORR.  With  taking  into  account  the  impact 
of  the  double  layer  capacitance  in  the  cathode  electrode,  the  overall 
dynamic  current  balance  is  thus  given  by: 


R2  =  fca20co  exp  (^-»7a) 


(23) 


where  0C o  represents  the  surface  coverage  of  platinum  catalyst 
with  adsorbed  CO  species.  Note  that  the  overpotential  rja  in  Eqs. 
(22)  and  (23)  is  the  absolute  difference  between  the  real  potential 
and  the  thermodynamic-equilibrium  potential  of  MOR  which  are 
referenced  to  a  standard  hydrogen  electrode  (SHE). 

The  overall  balance  of  0C o  is  then  given  by: 


d$co  _ 
dr  - 


-  ^2 


(24) 


Taking  into  account  the  current  generated  in  both  reaction  steps, 
the  kinetics  of  the  anodic  reductive  reaction  can  be  expressed  by 


Ja  =  (4fi3s)kal^i(l  -0co)exp  (°^r) a) 

Ci\/rQr*u  '  ' 


+{2FAs)ka2ecoexp(^^r1fj  (25) 

with  As  denoting  the  specific  surface  area  of  active  reaction  sites. 

On  the  DMFC  cathode,  the  kinetics  of  oxygen  reduction  reaction 
(ORR)  can  be  described  by  first-order  Tafel  expression 


where  Co2  represents  oxygen  concentration  in  the  gas  pores.  The 
term  kH,o2  denotes  the  Henry  factor  to  capture  the  effect  of  dissolv¬ 
ing  process.  The  overpotential  ijc  in  Eq.  (26)  represents  the  absolute 


Qcl  =  fceii  +  Wa  - 


jfccJcdxdy 

(Wc  +  Wr)/2 


(28) 


with  Jpara  denoting  the  ‘parasitic’  current  density,  which  can  be 
determined  by: 


gF  fW<+Wd/2  NulMEM/ccldx 
(Wc  +  Wr)/2 


(29) 


In  summary,  for  a  given  dynamic  cell  current  density,  the 
dynamic  anode  overpotential  can  be  determined  from  Eq.  (27), 
while  the  dynamic  cathode  overpotential  can  be  obtained  from  Eq. 
(28).  Finally,  the  dynamic  cell  voltage  can  be  obtained  from 


^Cell  =  Vo  -  ha 


he  -  kell 


^Contact  +^j 


(30) 


where  V0,  Rcontact  and  k  denote  the  thermodynamic  equilibrium 
voltage  of  a  DMFC,  the  contact  resistance  and  the  proton  conduc¬ 
tivity  of  the  membrane,  respectively. 


3.  Results  and  discussion 


The  above  described  governing  equations  for  the  cell  geometric 
dimensions  and  operating  parameters  listed  in  Table  2  subjected 
to  the  boundary  conditions,  along  with  the  constitutive  relations 
shown  in  Table  1  and  the  physicochemical  properties  listed  in 
Table  3,  are  solved  numerically  using  a  self-written  code,  which  was 
developed  based  on  the  SIMPLE  algorithm  with  the  finite-volume- 
method.  Note  that  the  present  transient  DMFC  model  is  extended 
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Table  2 

Cell  geometric  dimensions  and  operating  parameters 


Parameters 

Symbols 

Value 

Unit 

Anode  diffusion  layer  thickness 

^adl 

3.0  x  10-4 

m 

Anode  catalyst  layer  thickness 

5acl 

0.2  x  10-4 

m 

Membrane  thickness 

<5mem 

1.8  x  10-4 

m 

Cathode  diffusion  layer  thickness 

^cdl 

2.6  x  10-4 

m 

Cathode  catalyst  layer  thickness 

<5CC1 

0.2  x  10-4 

m 

Channel  width 

Wc 

1.0  x  lO"3 

m 

Rib  width 

Wr 

1.0  x  lO"3 

m 

Operation  temperature 

T 

333.15 

K 

Anode  inlet  pressure 

pin 

1.013  x  105 

Pa 

Cathode  inlet  pressure 

Pt 

1.013  x  105 

Pa 

Inlet  methanol  concentration  at  anode 

CM 

500 

mol  ITT3 

Inlet  methanol  vapor  concentration 

/-in 

UMV 

/-sat 

'“'MV 

mol  m3 

Inlet  oxygen  concentration  at  cathode 

c£ 

7.35 

mol  m3 

Inlet  water  vapor  concentration  at  anode 

CjSv 

/-sat 

WV 

mol  m3 

Inlet  water  vapor  concentration  at  cathode 

cjSv 

0 

mol  m3 

Inlet  liquid  saturation  at  anode 

sin 

1 

- 

Inlet  liquid  saturation  at  cathode 

Sin 

0 

- 

from  our  steady-state  DMFC  model  developed  earlier  [30,31  ],  which 
have  proved  to  be  competent  in  predicting  the  steady-state  mass 
transport  behaviors  in  a  liquid  DMFC. 


3. 1.  Dynamic  behavior  of  the  cell  voltage 

Fig.  2  shows  the  time-dependent  behavior  of  the  cell  voltage  in 
response  to  the  step  changes  in  current  density  from  100,  75  or 
50  mA cm-2  to  a  constant  value  of  10  mA cm-2.  As  seen  from  this 
figure,  in  response  to  each  change  in  current  density  the  cell  volt¬ 
age  sharply  rises  to  a  temporarily  higher  value  and  then  gradually 
lowers  to  a  stable  value.  Such  a  behavior  has  been  termed  as  the 
voltage  overshoot  resulting  from  a  dynamic  load  change.  Also,  it  is 
found  that  a  larger  step  change  in  the  cell  current  leads  to  stronger 
overshoot  behavior  in  the  cell  voltage.  With  respect  to  the  time 
characteristic  of  the  dynamic  voltage  response,  it  takes  over  30  s  for 
the  cell  voltage  to  reach  a  new  steady-state  for  all  the  cases  studied, 
closer  to  that  observed  experimentally  [23],  but  six  times  longer 
than  that  predicted  by  other  models  [23,24,27],  in  which  the  mass 
transport  processes  were  assumed  to  be  sufficiently  fast  such  that 
the  processes  can  be  treated  as  a  quasi-steady-state.  Apparently, 
this  assumption  of  fast  mass  transport  can  result  in  a  significant 
error  in  predicting  the  cell  dynamic  behavior.  The  time  charac¬ 
teristic  of  the  mass  transport  processes  will  be  discussed  later. 
Fig.  3  also  presents  the  voltage  overshoot  behaviors  in  response  to 
other  four  step  changes  in  current  density  from  60  to  10  mAcrrr2, 
from  70  to  20  mA cm-2,  from  80  to  30  mA cm-2,  and  from  90  to 


Table  3 

Physicochemical  properties 


Parameters 

Symbols 

Value 

Unit 

Ref. 

Porosity,  permeability 

ADL 

eadl-^adl 

0.7,  2.0  x  10-12 

-,  m2 

[33] 

ACL 

eacl  >  ^acl 

0.3, 1.0  x  10“14 

-,  m2 

MEM 

£mem>  kmem 

0.3,  2.0  x  10“18 

-,  m2 

CCL 

eccl  >  kcd 

0.3, 1.0  x  10“14 

-,  m2 

CDL 

ecdl  >  ^cdl 

0.7,  2.0  x  lO”12 

-,  m2 

[33] 

Nafion  volume  fraction  in  ACL 

eN,acl 

0.4 

_ 

MeOH  in  water 

^M.l 

10-5.4163-999.778 /T 

m2  s-1 

[14] 

MeOH  in  Nafion 

4.9  xlO-10  g[2436(i/333  —  l/r) 

m2  s-1 

[34] 

Methanol  vapor 

^M,g 

-6.954  X  10-6  +  4.5986 

m2  s_1 

[14] 

Diffusivities 

02  in  gas 

Do2,g 

x  10_8T+ 9.4979  x  10_11T2 

c  /  T  X  1.823 

1.775  x  10  5  (  275715  ) 

m2  s_1 

[14] 

Water  vapor  in  gas 

DWv,g 

2.56  x  10  (3077,5) 

m2  s-1 

[14] 

Viscosity  of  gas  phase 

Mg 

2.03  x  10-5 

kgm-1  s-1 

[36] 

Viscosity  of  liquid  phase 

Ml 

4.05  x  10-4 

kg  m_1s_1 

[35] 

Electro-osmotic  drag  coefficients  of  water  and 

nd,H20 

2.5 

- 

[37] 

methanol 

dd,M 

nd,H20*M 

- 

[37] 

Reaction  rate  constant  for  anode  reaction  step  1 

feat 

1.6  X  10-3 

mol  ITT2  s_1 

[1] 

Reaction  rate  constant  for  anode  reaction  step  1 

ka2 

OO 

b 

X 

0 

mol  m-2  s_1 

[1] 

Condensation  rate  constant  for  water 

kc 

5.0  x  10-5 

molatm-1  s-1  cm-3 

[13] 

Evaporation  rate  constant  for  water 

ke 

5.0  x  10-3 

atm^1  s-1 

[13] 

Henry  law  constant  for  oxygen 

ku,o2 

e(-666/T+ 14.1)  !RjT 

- 

[38] 

Henry  law  constant  for  methanol 

ku,M 

0.096eao4511(T-273) 

atm 

[14] 

Double  layer  capacitance  of  ACL 

Cacl 

1827 

Cm-2 

[23] 

Double  layer  capacitance  of  CCL 

CcCL 

907 

Cm-2 

[23] 

Surface  concentration  of  Pt  catalyst 

rpt 

0.11 

mol  m-2 

[1] 

Interfacial  transfer  rate  constant  for  methanol 

frig 

0.001 

m2  s_1 

[16] 

Specific  interfacial  area  between  liquid  and  gas 

Alg 

105 

m1 

[16] 

Specific  surface  area  of  anode  catalyst 

As 

5000 

m-1 

Proton  conductivity  in  membrane 

K 

7  3e[1268(l/298  -  1/T)] 

m-1 

[34] 

The  saturation  pressure  of  water  vapor 

lQgl0  P^sj 

-2.1794  +  0.02953 (T-  273)  -  9.1837 

atm 

[39] 

The  saturation  pressure  of  methanol  vapor 

nsat 

^MV 

x  10-5(T-  273 )3  + 1.4454  x  10-y(T-  273 )3 

atm 

[14] 

Thermodynamic  voltage 

Vo 

1.21 

V 

[14] 

Transfer  coefficient  of  anode  reaction  step  1 

«al 

0.5 

- 

[1] 

Transfer  coefficient  of  anode  reaction  step  2 

«a2 

0.5 

- 

[1] 

Transfer  coefficient  of  cathode 

ac 

1.0 

- 

[18] 

Cathode  exchange  current  density 

A  iref 

Av^o,o2 

1.14  x  103 

A  ITT3 

[18] 

Anode  reference  concentration 

CMf 

1000 

mol  ITT3 

[1] 

Cathode  reference  concentration 

r  ref 

lo2 

0.52 

mol  m-3 

[18] 
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Time  /  s 

Fig.  2.  Dynamic  behavior  of  the  cell  voltage  in  response  to  the  current  density 
change  from  (a)  100  to  10  mA cm-2,  (b)  75  to  10mAcm_2,and(c)50to  lOmAcirr2. 

40mA cm-2.  Note  that  in  all  the  cases,  the  decrement  in  current 
density  is  kept  the  same  as  50  mA  cm-2.  Interestingly,  the  voltage 
overshoot  becomes  weaker  with  increasing  the  initial  current  den¬ 
sity.  This  result  is  qualitatively  consistent  with  the  experimental 
findings  [23]. 

In  order  to  clarify  the  dominant  factor  that  causes  the  voltage 
overshoot  in  response  to  a  change  in  current  density,  the  time- 
dependent  behaviors  of  both  the  anode  and  cathode  overpotentials 
are  plotted  in  Fig.  4.  Clearly,  in  response  to  the  decrease  in  cur¬ 
rent  density,  the  anode  overpotential  drops  suddenly  from  a  higher 
level  to  a  lower  level  for  all  the  cases  studied,  due  primarily  to  the 
decrease  in  the  charge  transfer  resistance  of  the  anode  electrode 
reaction.  Unlike  the  anode  overpotential,  the  cathode  overpoten¬ 
tial  shows  a  significant  undershooting,  which  is  believed  to  be  the 
main  reason  that  results  in  the  cell  voltage  overshoot.  Also,  it  can  be 
found  that  it  takes  about  30  s  for  the  anode  overpotential  to  reach 
a  new  steady-state,  whereas  it  takes  almost  60  s  for  the  cathode 
overpotential  to  reach  a  new  steady-state  after  the  current  den¬ 
sity  changes.  The  different  time  characteristics  imply  the  different 
physicochemical  processes  occurring  in  each  electrode.  In  the  fol¬ 
lowing,  we  discuss  about  the  effects  of  different  physicochemical 
processes  on  the  cell  dynamic  operation  characteristics. 


Time  /  s 

Fig.  3.  Dynamic  behavior  of  the  cell  voltage  in  response  to  the  current  density 
change  of  (a)  60  to  10  mAcirr2,  (b)  70  to  20  mAcirr2,  (c)  80  to  30  mAcirr2,  and  (d) 
90  to  40  mA  cm-2. 


Fig.  4.  Dynamic  behaviors  of  (a)  anode  overpotential  and  (b)  cathode  overpotential 
in  response  to  the  change  in  current  density. 

3.2.  The  impact  of  the  double  layer  capacitance 

To  make  clear  the  influence  of  the  double  layer  capacitance,  we 
compare  (a)  the  model  that  takes  account  the  impact  of  the  DLC 
and  (b)  the  model  that  ignores  the  impact  of  the  DLC.  The  results 
are  shown  in  Fig.  5.  Clearly,  the  voltage  overshoot  becomes  weaker 
with  the  consideration  of  DLC.  This  result  suggests  that  the  DLC  can 
smoothen  the  change  in  voltage  in  response  to  the  change  in  cell 
current.  Also,  it  is  noticed  that  the  impact  of  the  DLC  is  transitory, 
which  only  reacts  on  the  dynamic  voltage  within  a  few  seconds  after 
the  current  density  changes  because  the  DLC  charging/discharging 
process  is  generally  very  fast. 

3.3.  Surface  coverage  of  CO  on  the  Pt  catalyst  in  the  ACL 

As  indicated  in  Eq.  (25),  the  adsorbed  CO  species  on  the  Pt  cat¬ 
alyst,  0co,  known  as  catalyst  poisoning,  is  closely  related  to  the 
kinetics  of  MOR  and  thus  may  play  a  role  on  the  anode  performance 
as  the  current  density  changes.  Fig.  6  shows  the  time-dependent 
variation  of  0C o  in  response  to  the  changes  in  current  density.  Nor¬ 
mally,  0Co  is  relatively  small  when  the  fuel  cell  discharges  at  a 
stably  high  current  density.  As  the  current  density  suddenly  drops, 
there  are  more  free  Pt  catalyst  sites  than  needed  for  the  complete 
methanol  oxidation,  so  methanol  adsorption  and  the  dehydro¬ 
genation  step  become  prevailing,  which  consequently  results  in  an 
increase  in  0CO.  Clearly,  as  seen  from  Fig.  6,  the  surface  coverage 
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Fig.  5.  Dynamic  behaviors  of  the  anode  overpotential,  the  cathode  overpotential  and 
the  cell  voltage  in  response  to  the  change  in  current  density  from  100  to  10  mA  cm-2 . 
Case  A:  with  the  double  layer  capacitance;  Case  B:  without  the  double  layer  capac¬ 
itance. 


of  CO  species  on  Pt  catalyst  increases  from  around  0.44  to  around 
0.84  when  the  current  density  drops  from  100  to  10  mA  cm-2.  With 
respect  to  the  time  characteristic  of  the  catalyst  poisoning  process, 
it  takes  over  60  s  for  0Co  to  reach  a  new  steady-state.  However,  it 
is  noticed  that  although  there  has  been  large  variation  in  0C o  after 
the  current  density  changes,  the  anode  overpotential  exhibits  a  rel¬ 
atively  constant  value  for  the  case  without  consideration  of  double 
layer  capacitance,  as  shown  in  Fig.  5.  This  indicates  that  the  anode 
overpotential  is  relatively  insensitive  to  the  change  in  surface  cov¬ 
erage  of  CO  species  on  Pt  catalyst  in  the  ACL.  This  result  is  consistent 
with  the  work  by  Krewer  and  Sundmacher  [23]. 

3.4.  Mass  transport  of  methanol  in  the  DMFC 

Another  important  parameter  that  may  influence  the  kinetics 
of  the  anode  electrode  reaction  is  the  methanol  concentration  in 
the  ACL.  Fig.  7  shows  the  time-dependent  variations  in  methanol 
concentration  and  gas  void  fraction  in  the  ACL  in  response  to  the 
changes  in  current  density.  As  the  current  density  is  suddenly 
decreased,  the  demand  of  the  methanol  for  electro-oxidation  is 
decreased,  thus  breaking  the  original  balance  of  methanol  in  the 


Fig.  6.  Dynamic  behavior  of  the  CO  surface  coverage  (0co)  in  the  ACL  in  response  to 
the  change  in  current  density. 


Fig.  7.  Dynamic  behaviors  of  the  methanol  concentration  and  gas  void  fraction  in 
the  ACL  in  response  to  the  change  in  current  density. 


anode  electrode.  Consequently,  the  methanol  begins  to  accumu¬ 
late  in  the  ACL  until  it  reaches  a  new  steady-state,  as  can  be  seen 
in  the  Fig.  7.  Finally,  the  methanol  concentration  in  the  ACL  reaches 
the  same  steady-state  value  because  the  final  current  densities  for 
all  cases  studied  are  the  same.  Generally,  it  takes  over  90  s  for  the 
methanol  concentration  in  the  ACL  to  reach  a  new  steady-state, 
showing  that  the  mass  transport  of  methanol  is  sluggish  in  the  liq¬ 
uid  feed  DMFC.  It  should  be  pointed  out  that  although  the  methanol 
concentration  in  the  ACL  is  greatly  changed  after  the  current  den¬ 
sity  drops,  it  does  not  show  a  significant  influence  on  the  anode 
overpotential,  probably  because  the  kinetics  of  the  anodic  MOR 
is  insensitive  to  methanol  concentration  once  it  exceeds  a  certain 
value  [14,32].  Similar  to  the  methanol  transport,  the  transport  of 
the  gas  in  the  ACL  is  also  relatively  slow.  As  shown  in  the  figure, 
it  takes  approximate  60s  for  the  gas  void  faction  in  the  ACL  to 
reach  a  new  steady-state.  To  have  a  deeper  understanding  of  the 
dynamic  characteristics  of  the  methanol  transport  in  the  DMFC, 
the  time-dependent  variations  in  methanol  concentration  in  dif¬ 
ferent  regions  inside  the  cell  are  displayed  in  Fig.  8.  The  methanol 
concentrations  at  six  representative  regions  inside  the  DMFC  are 
probed,  i.e.:  the  middle  points  of  the  ADL,  the  ACL  and  the  mem¬ 
brane  under  the  midlines  of  the  channel  and  the  rib.  As  can  be  seen 
from  Fig.  8,  the  methanol  concentration  in  each  region  increases  in 
response  to  the  sudden  drop  in  current  density.  In  particular,  the 
variation  in  methanol  concentration  in  the  ACL  is  much  steeper  and 
larger  than  in  the  other  regions.  Generally,  it  takes  about  45  s  for  the 
methanol  concentration  under  the  midline  of  the  channel  to  reach 
a  new  steady-state,  as  shown  in  Fig.  8a.  In  contrast,  almost  3  min 
is  needed  for  the  methanol  concentration  under  the  midline  of  the 
rib  to  reach  a  new  steady-state,  as  shown  in  Fig.  8b.  The  longer 
response  of  the  methanol  concentration  in  the  regions  under  the 
rib  is  resulted  from  the  large  transfer  resistance  and  longer  trans¬ 
fer  distance  for  the  methanol  in  those  regions.  The  evolution  of 
the  methanol  concentration  profile  across  the  membrane  electrode 
assembly  under  the  centerline  of  the  channel  is  displayed  in  Fig.  9. 
As  can  be  seen  from  this  figure,  the  methanol  concentration  in  the 
ACL  and  its  adjacent  regions  responds  very  fast  within  the  first  sec¬ 
ond  after  the  current  density  drops,  whereas  in  other  regions  far 
away  from  the  ACL  the  methanol  concentration  remains  almost 
the  same.  This  hysteresis  is  attributed  to  the  slow  diffusion  of  the 
methanol  through  the  ADL  and  the  MEM. 

The  variations  in  the  methanol  concentration  in  the  ACL  and 
in  the  MEM  can  significantly  influence  the  flux  of  methanol  per¬ 
meation  through  the  membrane  to  the  cathode.  Fig.  10  shows  the 
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Fig.  8.  Dynamic  behavior  of  the  methanol  concentration  in  anode  electrode  and 
membrane  in  response  to  the  change  in  current  density  from  100  to  10  mAcm-2. 

time-dependent  behavior  of  the  methanol  crossover  in  response 
to  current  density  changes.  Note  that  flux  of  methanol  crossover  is 
expressed  in  terms  of  the  ‘parasitic’  current  density.  As  can  be  seen, 
the  flux  of  methanol  crossover  gradually  increases  in  response  to 
the  drop  in  current  density.  It  takes  over  90  s  for  the  flux  of  methanol 


Fig.  9.  Dynamic  evolution  of  the  methanol  concentration  profiles  through  the  MEA 
under  the  centerline  of  the  channel  in  response  to  the  change  in  current  density 
from  100  to  10 mAcm-2. 


Fig.  10.  Dynamic  behavior  of  the  methanol  crossover  in  response  to  the  change  in 
current  density. 


crossover  to  reach  a  new  steady-state.  The  time  characteristics  of 
the  methanol  crossover  are  obviously  in  the  same  time  range  as  the 
methanol  concentration  increase  in  the  ACL.  As  concluded  earlier, 
the  dynamic  methanol  crossover  results  in  the  dynamic  response  of 
cathode  mixed  overpotential,  which  has  proved  to  be  the  dominant 
reason  that  results  in  the  cell  voltage  overshoot  behavior. 

3.5.  Mass  transport  of  oxygen  and  water  in  the  CCL 

Both  the  oxygen  concentration  and  water  saturation  in  the  CCL 
change  with  the  current  density,  which  in  turn  may  affect  the  cath¬ 
ode  performance.  Fig.  11  shows  the  time-dependent  variations  in 
oxygen  concentration  and  liquid  saturation  in  the  CCL  in  response  to 
the  changes  in  current  density.  As  can  be  seen,  the  oxygen  concen¬ 
tration  increases  to  a  relatively  higher  value  within  a  few  seconds 
resulting  from  the  sudden  decrease  in  the  current  density.  Then,  it 
gradually  decreases  and  reaches  a  steady-state.  The  whole  process 
takes  over  90  s.  With  respect  to  the  liquid  saturation  in  the  CCL,  it 
gradually  decreases  to  a  new  steady-state  in  response  to  the  cur¬ 
rent  density  drop.  Although  both  the  oxygen  concentration  and  the 
liquid  saturation  in  the  CCL  change  with  a  sudden  drop  in  current 
density,  the  variation  is  relatively  small  for  operating  conditions 
considered.  However,  it  should  be  mentioned  that  the  liquid  sat¬ 
uration  in  the  CCL  may  have  a  significant  impact  on  the  cathode 
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Fig.  11.  Dynamic  behaviors  of  the  oxygen  concentration  and  liquid  water  saturation 
in  the  CCL  in  response  to  the  change  in  current  density. 
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performance  wherever  the  water  flooding  is  becoming  serious  in 
the  DMFC  cathode  electrode. 

4.  Conclusions 

A  transient,  two-dimensional,  two-phase  mass  transport  model 
is  developed  to  investigate  the  dynamic  operating  behavior  of  a  liq¬ 
uid  feed  DMFC.  Various  processes  that  may  affect  the  cell  operating 
dynamics  are  studied.  The  numerical  results  show  that  the  transient 
voltage  exhibits  a  significant  overshoot  in  response  to  a  sudden 
change  in  current.  The  degree  of  the  overshoot  depends  on  the  mag¬ 
nitude  of  the  change  in  current  density  and  the  value  of  the  initial 
current  density.  It  is  found  that  the  methanol  permeation  through 
the  membrane  to  the  cathode  results  in  a  strong  cathode  overpo¬ 
tential  overshoot,  which  is  the  predominated  reason  that  leads  to 
the  voltage  overshoot  behavior.  In  contrast,  the  anode  overpotential 
is  found  to  be  insensitive  to  the  change  in  methanol  concentration 
and  CO  surface  coverage  in  the  ACL.  The  effect  of  double  layer  capac¬ 
itance  on  the  cell  dynamic  operation  is  also  studied  and  the  results 
show  that  the  DLC  can  smoothen  the  change  in  voltage  in  response 
to  the  change  in  cell  current.  The  dynamic  response  of  mass  trans¬ 
port  to  a  change  in  current  density  is  found  to  be  rather  slow.  In 
particular,  it  is  shown  that  the  slow  mass  transport  of  methanol  is 
one  of  the  key  factors  that  influence  the  cell  dynamic  operation. 
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